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To understand effects of orbital degeneracy on magnetism, in particular effects of Hund's rule 
coupling, we study the two-orbital Hubbard model on a square lattice by a variational Monte Carlo 
method. As a variational wave function, we consider a Gutzwiller projected wave function for a 
staggered spin and/or orbital ordered state. We find a ferromagnetic phase with staggered orbital 
order around quarter-filling, i.e., electron number n — 1 per site, and an antiferromagnetic phase 
without orbital order around half-filling n — 2. In addition, we find that another ferromagnetic 
phase without orbital order realizes in a wide filling region for large Hund's rule coupling. These 
two ferromagnetic states are metallic except for quarter filling. We show that orbital degeneracy 
and strong correlation effects stabilize the ferromagnetic states. 

PACS numbers: 71.10.Fd, 75.30.Kz 



Mechanism of itinerant ferromagnetism is a long stand- 
ing problem in physics of condensed matter. As a sim- 
ple model for itinerant ferromagnetism, the single-orbital 
Hubbard model has been studied intensively, but it has 
been revealed that it is difficult to stabilize a ferromag- 
netic state in the Hubbard model with only nearest- 
neighbor hopping on simple lattices such as a square lat- 
tice. 

One possible improvement to the Hubbard model for 
ferromagnetism is a modification of the band structure. 
Since the early stage of the study of ferromagnetism, 
a large density of states around Fermi level, as in an 
fee lattice with appropriate filling, has been suggested 
to stabilize a ferromagnetic stateji' 2 -' 3 - Indeed, realization 
of ferromagnetic ground states is proven for some flat 
band systems^ and nearly flat band systems^ It is also 
shown that ferromagnetism occurs for an fec-type infi- 
nite dimensional lattice and for an fee lattice by using a 
dynamical mean-field theory, - while a dynamical mean- 
field theory for a hypercubic lattice does not show fer- 
romagnetism^ For finite dimensions, it has been shown 
that ferromagnetism can occur by including next-nearest 
hopping, which induces Van Hove singularity, for a chain 9 
and for a square latticed 

Another possible improvement is inclusion of orbital 
degree of freedom, which may be important to deal with 
realistic situations in transition metals. For orbitally de- 
generate systems, it has been suggested that intra-atomic 
Hund's rule coupling can stabilize f er romagnet ism i 1 1 ' 1 2 1 1 3 
The simplest extended model of the single-orbital Hub- 
bard model including orbital degree of freedom is the 
two-orbital Hubbard model. This model shows ferromag- 
netism with antiferro-orbital order at quarter-filling, i.e., 
electron number n = 1 per site, in the strong Coulomb 
interaction limits 

This ferromagnetic state at n = 1 is insulating. Thus, 
it is an interesting problem as to what extent is the ferro- 
magnetic state stable against doping of electron or hole 
which makes the system metallic. It is found that the 



ferromagnetic state is stable to some extent against dop- 
ing in one dimension^ and in infinite dimensions i 16 > 17 
For other finite dimensions, there are few studies on 
doping effects on magnetism of the two-orbital Hubbard 
model. Sakai et al^ have studied the two-orbital Hub- 
bard model on an fee lattice by a dynamical mean-field 
theory, and have stressed importance of the lattice struc- 
ture and Hund's rule coupling for ferromagnetism. 



To understand magnetism of the two-orbital Hubbard 
model deeply, we have to investigate the model for many 
parameter sets, since the two-orbital Hubbard model has 
a parameter for Hund's rule coupling in addition to that 
for the Coulomb interaction. However, such an extensive 
study has been difficult for a two-orbital model beyond 
the Hartree-Fock approximation, since there are as many 
as sixteen electron configurations at each site and it is 
hard for numerical calculations. 



In this paper, to overcome such difficulty, we apply 
a variational Monte Carlo methodic to the two-orbital 
Hubbard model on a square lattice. We use a Gutzwiller 
projected wave function as a variational wave function. 
This wave function is simple enough but includes correla- 
tion effects, and we can evaluate energy for several states 
and for various parameters. In particular, we can con- 
struct a phase diagram by varying the value of Hund's 
rule coupling and filling n. Thus, we can investigate over- 
all feature of the two-orbital Hubbard model. At n = 1 , a 
similar model without considering orbital order— and the 
two-orbital Hubbard model considering possibility of or- 
bital order— have been studied by the variational Monte 
Carlo method, but the effect of doping has not been in- 
vestigated by these studies, which is a main topic of the 
present paper. 
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The two-orbital Hubbard model is given by 
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where Ci T(T is the annihilation operator of the electron 
at site i with orbital t (= 1 or 2) and spin a (=t or 
I), CkTCT is the Fourier transform of it, n lTa — cj TCT c JT(T , 
and rii T = ^2,(y n 'ira- The coupling constants U, U', J, 
and J' denote the intra-orbital Coulomb, inter-orbital 
Coulomb, exchange, and pair-hopping interactions, re- 
spectively. We use the relation U = U' + J + J' , which 
is satisfied in several orbital-degenerate models such as 
a model for p-orbitals, a model for e g orbitals, and a 
model for t2 g orbitals.— We also use the relation J = J', 
which holds if we can choose wave functions of orbitals 
reali^ We consider only a nearest-neighbor hopping in- 
tegral t for both orbitals, and the kinetic energy is given 
by £k = 2i(cosfca; + cos ky). Here we have set the lattice 
constant unity. 

We consider the variational wave function given by 
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where Pq is the Gutzwiller projection operator for on- 



site density correlation 
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denotes projec- 



tion onto the state 7 at site i and g 1 is the variational 
parameter controlling the probability of state 7. There 
are sixteen states at each site in the present two-orbital 
model. The Hartree-Fock type wave function |$), which 
describes a charge, spin, orbital, and spin-orbital coupled 
ordered state, is given by 



i*>= n 



|0>, 
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where |0) is the vacuum. The quasiparticles occupy N a 
states for each spin a from the lowest quasiparticle energy 
state, where N a is the number of electrons with spin a. 
The energy of the quasiparticle in the ordered state is 
given by 



,0) 
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The creation operators of quasiparticles are given by 

b Va = U ^raci T<J + Sgn(A T<7 )w kT a4+Qr<7' ( 5 ) 

where Q = (it, it) is the ordering vector considered in this 
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The quasiparticle gap in the ordered state is given by 
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where A c , A 
spin, orbital 



A , and A so denote the gaps for charge, 
and spin-orbital ordered states, respec- 
tively, and we take them as variational parameters. 

Here, we have chosen z-component of spin for the or- 
dered state. We can choose x- or y-component, but they 
are equivalent due to the rotational symmetry in the spin 
space. On the other hand, there is rotational symmetry 
only in the z-x plane in the orbital space. For orbital 
order, in addition to z-component as in Eq. ([9]), we have 
also investigated possibility of y-component order. The 
model Hamiltonian |T]) can be rewritten in terms of ba- 
sis states of y-component of orbital by replacing interac- 
tion parameters with tilde: U = (U + V + J — J')/2, 
U' = (U + V - J + J')/2, J = (U -U' + J+ J')/2, 
and J' = (-Z7 + U' + J + J')/2. Thus we can study 
y-component orbital order with the same form for the 
variational function by simply changing interaction pa- 
rameters. 

We evaluate the expectation value of energy for the 
variational wave function by using the Monte Carlo 
method, and optimize these gap parameters and the 
Gutzwiller parameters to minimize energy. For the opti- 
mization we use a fixed sampling method i 26 ' 27 The num- 
ber of parameters in the Gutzwiller projection can be 
reduced from sixteen to ten by considering equivalence 
of the two orbitals 3^ We can further reduce the number 
of the Gutzwiller parameters to seven when we consider 
conservation of spin, i.e., when we fix N a . Thus the to- 
tal number of the variational parameters is eleven. Note 
that, in this study, we take partial ferromagnetic states 
into consideration and spin j and J, are not equivalent 
for these states. Thus we cannot reduce the number of 
the parameters further by considering spin states. We 
can also evaluate energy by fixing some parameters, for 
example, we set all the gap parameters zero for a para- 
magnetic state. The calculations have been done for a 
12 x 12 lattice with periodic- antiperiodic boundary con- 
ditions. 

Figure [1] shows energy E in ordered states per site 
measured from energy -Ep ar a per site in the paramag- 
netic state as functions of n for U/t = 15 and J/t = 2 
as an example. Statistical errors are much smaller than 
the symbol sizes. In the low filling region, the param- 
agnetic state is stable. Around quarter-filling n = 1, 
the ferromagnetic phase with antiferro-orbital order ap- 
pears. In the filling region 1.38 < n < 1.78, the ferro- 
magnetic phase without orbital order appears. Around 
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FIG. 1: (Color online) Filling dependence of energy for several 
states measured from that of the paramagnetic state: Spin- 
AF Orbital- Para (m = 0, A s ^ 0, A D = 0, solid triangles), 
Partial-Ferro m = 0.5 Orbital-Para (A s = 0, A = 0, solid 
squares), Spin-Ferro Orbital-Para (m = 1, A s = 0, A = 0, 
solid circles), and Spin-Ferro Orbital-AF (m = 1, A s = 0, 
A 7^ 0, open circles). 



half- filling n = 2, antifcrromagnetic phase appears due to 
nesting of the Fermi surface. While we have calculated 
energy of partial ferromagnetic states, for example, for 
m = (JVf - N l )/(N 1 + N L ) = 0.5 shown in Fig. □ these 
states do not become the ground state. 

Here we comment on the ground states at n — 1 and 
2. The chemical potential /x can be obtained from /i = 
dE/dn in the ground state and we also obtain the relation 



dn d 2 E 
d^ji dn 2 



= 1. 



(10) 



Thus, in an insulating state, i.e., dn/d/i = 0, the second 
derivative of energy E with respect to n should diverge 
and vice versa. At n — 1, the ground state energy has 
a cusp as shown in Fig. [1] and the ground state is insu- 
lating. Note that we can obtain energy for n > 2 from 
the present data by using electron-hole symmetry of the 
model, and we find a cusp in the ground state energy also 
at n = 2. 

Figure f2] shows phase diagrams for U/t = 9 and 
15. First, we explain the phase diagram Fig. [UJb) for 
U/t = 15. The ferromagnetic phase with antiferro-orbital 
order appears around quarter-filling as is expected from 
the effective Hamiltonian for the strong coupling limit. 
However, it is found that this ferromagnetic phase is eas- 
ily destabilized by doping. To stabilize this ferromag- 
netic phase in a wider filling region, much larger value of 
Coulomb interaction is necessary. 

This ferromagnetic phase is also destabilized by in- 
creasing Hund's rule coupling J. Hund's rule coupling is 
usually expected to stabilize magnetically ordered states, 
but the effective interaction U c s = U' — J between dif- 
ferent orbitals is reduced by Hund's rule coupling. Thus, 



0.3 - 

0.2 - 

0.1 - 




U/t=9 Spin-Ferro 
Orbital-Para 



Para 




Spin-AF 
4 Orbital-Para 



(a) 



0.3 - 
0.2 - 



U/t=\5 



Spin-Ferro 
vOrbital-Para 



Para 



Spin-Ferro 
0.1 h Orbital-AF"" 



(b) 




<Spin-AF 
\Orbital-Para 



0.5 



1 1.5 



FIG. 2: (Color online) Phase diagrams for (a) U/t = 9 and 
for (b) U/t — 15. Solid lines denote first-order transitions and 
dashed lines denote second-order transitions. 



Hund's rule coupling destabilizes the orbital order, and as 
a result, the ferromagnetic state supported by the orbital 
order is also destabilized. Note that in the ferromagnetic 
phase, there is rotational symmetry in the orbital space, 
and orbital order for x-, y-, and z-components are equiv- 
alent. 

At higher filling region, another ferromagnetic phase 
without orbital order appears in a large Hund's rule cou- 
pling region. The ferromagnetic phase extends in a wide 
parameter region for n > 1, while not for n < 1. This 
finding is in agreement with the statement that double- 
exchange type mechanism works well for n > 1 since the 
probability of double occupancy is high, but it is less 
effective for n < 1.— 

In this ferromagnetic phase at n ~ 1 and J/U ~ 0.3, 
we expect orbital-antiferro order, since in the ferromag- 
netic state the model is reduced to the single-orbital Hub- 
bard model with effective interaction U e R if we regard 
spin in the single-orbital Hubbard model as orbital and 
around n = 1 an orbital-antiferro state should occur. 
However, it is difficult to distinguish a small energy dif- 
ference between orbital-para and orbital-antiferro states 
around there due to a small value of U c g. Thus, the phase 
boundary between spin-ferro orbital-para and spin-ferro 
orbital-antiferro in Fig. EJb) is merely a eye guide. 

Around half-filling, the antiferromagnetic phase ap- 
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pears as is expected from the nesting of Fermi surface. 
The phase transition from the paramagnetic phase to 
the antiferromagnetic phase is second-order. We have 
checked that the energy difference between these phases 
is proportional to (n — n c ) 2 for n > n c , where n c is the 
critical filling. 

Note that the spin-ferro orbital-antiferro state at n = 1 
and the spin-antiferro orbital-para state at n = 2 are 
insulating, and other ground states are metallic. 

By reducing the Coulomb interaction U, the regions of 
the ordered phases become narrower as shown in Fig.[2Ja) 
for U/t = 9. In particular the spin-ferro orbital-antiferro 
state disappeared. At n = 1, the ferromagnetic state 
with orbital order disappears at U/t < 10.— The other 
ferromagnetic phase without orbital order is also eas- 
ily destabilized by reducing the Coulomb interaction. 
This fact indicates that realization of ferromagnetism is a 
strong correlation effect. On the other hand, the antifer- 
romagnetic phase around half-filling, which is stabilized 
by the nesting of the Fermi surface and can be obtained 
with a weak-coupling theory, realizes in a wide region 
even at U/t — 9. 

Note that we have also calculated energy of the single- 
orbital Hubbard model within the Gutzwiller wave func- 
tion, and we have found that a much larger value of 
Coulomb interaction U/t > 23 is necessary to stabilize 
a ferromagnetic phase. Thus, the orbital degeneracy and 



Hund's rule coupling are important ingredients for real- 
ization of ferromagnetism with a moderate value of the 
Coulomb interaction. 

To summarize, we have studied the two-orbital Hub- 
bard model on a square lattice by a variational Monte 
Carlo method. We have considered charge, spin, orbital, 
and spin-orbital coupled ordered states for the variational 
wave function. Then, we have constructed phase dia- 
grams for the ground states. We find a narrow region 
of the ferromagnetic state with orbital-antiferro order 
around quarter-filling and a wide region of ferromagnetic 
phase without orbital order at large Hund's rule coupling 
for U/t = 15. The ferromagnetic phase with orbital or- 
der is easily destabilized by doping and by reducing the 
Coulomb interaction. The ferromagnetic phase without 
orbital order is also destabilized strongly by reducing the 
Coulomb interaction. Thus, realization of ferromagnetic 
states is a strong correlation effect. Investigation of ef- 
fects of realistic anisotropic hopping integral depending 
on orbital and further improvement of the variational 
wave function arc important future problems. 
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